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ABSTRACT 

Obscured or narrow-line active galaxies offer an unobstructed view of the quasar environment in the presence 
of a luminous and vigorously accreting black hole. We exploit the large new sample of optically selected lumi- 
nous narrow-line active galaxies from the Sloan Digital Sky Survey at redshifts 0.1 < z < 0.45, in conjunction 
with follow-up observations with the Low Dispersion Survey Spectrograph (LDSS3) at Magellan, to study the 
distributions of black hole mass and host galaxy properties in these extreme objects. We find a narrow range in 
black hole mass ((log M^yiIMq) = 8.0 ±0.7) and Eddington ratio ((log Lboi/^Edd) = —0.7 ±0.7) for the sample as a 
whole, surprisingly similar to comparable broad-line systems. In contrast, we infer a wide range in star formation 
properties and host morphologies for the sample, from disk-dominated to elliptical galaxies. Nearly one-quarter 
have highly disturbed morphologies indicative of ongoing mergers. Unlike the black holes, which are apparently 
experiencing significant growth, the galaxies appear to have formed the bulk of their stars at a previous epoch. On 
the other hand, it is clear from the lack of correlation between gaseous and stellar velocity dispersions in these 
systems that the host galaxy interstellar medium is far from being in virial equilibrium with the stars. While our 
findings cast strong doubt on the reliability of substituting gas for stellar dispersions in high luminosity active 
galaxies, they do provide direct evidence that luminous accreting black holes influence their surroundings on a 
galaxy-wide scale. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert 
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1. introduction 

Over recent years the cosmological status of accreting super- 
massive black holes (BHs) has been elevated from an extreme 
rarity to a ubiquitous phase in the life cy cle of all massive galax- 
ies (e.g., Kormendv & Richstone 1995). There appears to be di- 
rect evidence for a connection between BH and galaxy growth 
in the tight scaling relations between BH mass and bulge 
properties (e.g., the Mbh-c* relation; Tremai ne et all 12002). 
Nonetheless, all of the classic conundrums regarding BH re- 
main. There is not yet a reliabl e measurement of the length 
of typical BH growth episodes (Martini 200 4!), or the typi- 
cal mode of BH growth (e.g.jKollmeier et al.l l2006; Shen et al.l 
2008| iGavignaud et alj|2008l) . Furth ermore, while it has been 



argued both on obse rvatio nal (e.g., |S anders & Mirabel 1996; 
ICanalizo & Stocktonll2000l: iGonzalez Delgado et al.ll200ll) and 
theoretical dMihos & Hernquistll994l;lGoodmanl2003h grounds 
that star formation and accretion activity are temporally coin- 
cident, quantitative stellar population constraints for large sam- 
ples of accreting B Hs at high luminosities are notoriously dif- 
ficult to obtain (e.g.jBoroson & Okdll984j; | Bahcall et alJll997t 
Cana lizo & Stockton! 1200 it iLetawe et al.l 120071 ; iJahnke et al.l 
2007). The accretion power from a luminous active galactic nu- 



cleus (AGN) may easily outshine the surrounding host galaxy 
by a factor of 100, which makes studies of AGN host galaxy 
stellar populations and morphologies a considerable challenge. 



Occasionally nature is particularly kind and puts obscur- 
ing material between us and the accretion disk. In these 
cases, we know that an accreting BH is present due to the 
ratios of strong emissi on lines that impl y a very hard ion- 
izing spectrum (e.g.. iBaldwin etaTl [l981). and BHs selected 
this way are known as narrow-line or Type 2 AGNs. These 
objects are an ideal tool with which to study AGN host 
galaxies, since the galaxy to AGN contrast ratio is maxi- 
mized. Much of what we know about the stellar popula- 
tions of AGN host galaxies is derived from narrow-line ob- 
jects (e.gjHeckman et al.lll997l;lGonzalez Delgado et al.ll200lh 
Cid Fernande s et alj 120041) . However, until recently there 
were not large samples of homogeneously selected, luminous 
narrow-line AGNs. The situation has been remedied by the 
large spectrosc opic database pro v ided by the Sloan Digital Sky 
Survey (SDSS: lYork etai1l2000h . IZakamska et ail (l2003h used 
the strong and ubiquitous [O III] A5007 line emission as a proxy 
for bolometric luminosity in order to select a sample of narrow- 
line AGNs as luminous as any known broad-li ne active galaxies 
with z < L More recently Reves et al. (2008) presented a sam- 
ple of nearly 1000 narrow-line AGNs with z < 0.8 and bolo- 
metric luminosities as high as 10 47 erg s" 1 or more. 

We have initiated a program to obtain deeper follow-up long- 
slit spectroscopy of the Reyes et al. sample at z < 0.45. Our pri- 
mary goal is to study the evolutionary state of the host galaxies 
of vigorously accreting BHs. Both stellar velocity dispersions 
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(a») and stellar population age can be measured from the stel- 
lar continua. From 0% we derive not only a mass-scale for the 
galaxy, but also an estimate of the BH mass and in turn Edding- 
ton fractions for the sample, assuming that theMBH-c* relation 
holds (which may not be the case for active ga laxies at these 
redshifts; e.g.. lWoo et aT]l2006t iKim et alj F2008). By compar- 
ing the structural properties of the hosts to those of inactive 
galaxies at a similar mass, combined with our measurements of 
stellar ages, we can hope to determine whether the host galaxies 
are evolving coevally with their BHs. Of course, as we discuss 
in §7, it is not yet clear whether broad- and narrow-line AGNs 
are truly identical populati ons seen at differ ent orientations (as 
suggested by unification; lAntonuccl 1 19931) . or whether they 
represent an evolutionary sequence (as suggested by the dif- 
ferences in star formation properties , among other things, e .g., 
IKim et alJl2006tlHopkins etatll2006l:IZakamska et al.ll2008l) . 

Throughout we assume the following cosmological param- 
eters to calculate distances: Ho = 100 h = 70 km s" 1 Mpc -1 , 
n m = 0.30, and n A = 0.70. 

2. THE SAMPLE, OBSERVATIONS, AND DATA REDUCTION 

Our goal was to select nearby targets, for which we could 
hope to resolve the emission-line structure, that are neverthe- 
less at the upper end of the luminosity function. We thus lim- 
ited our attention to targets with z < 0.45, such that [O III] 
is accessible in the observing windows (see below) and with 
L[om] > 10 42 erg s" 1 . While arbitrary, our luminosity cut is cho- 
sen to match the classic definition of a quasar. Altogether, we 
obtained spectroscopic observations for 15 targets. Our sample 
spans, by construction, a much narrower range in redshift (and 
thus luminosity) than the parent sample. 

All of the data discussed in this paper were collected over 
two observing runs using the Low Dispersion Survey Spec- 
trograph (LDSS3) at the Clay-Magellan telescope on Las 
Campanas. LDSS3 is a wide-field (8.3' diameter) imaging 
spectrograph, which was up graded from the original LDSS2 
dAHington-Smifh et al.l fl994) to have increased sensitivity in 
the red. For each target we obtained 2-5 min acquisition im- 
ages in g and r, and sometimes i. We then spent at least 1 hr per 
target at the primary slit position, as well as additional time as 
noted in Table 1, using a l"x4' slit. Slit position angles were 
chosen by manual inspection of the SDSS color images, with 
an eye to investigate color gradients or low-surface brightness 
material (Table 1). The seeing ranged from 0"6 to 1 "5 over the 
two runs, but was typically 1"; with a plate scale of 0."189 the 
spectra are well sampled in the spatial direction. Two spectro- 
scopic settings were required in order to ensure that [O III] was 
observed over the full redshift range of our sample. The pri- 
mary setting was the reddest position of the VPH-Blue grism, 
which covers 4300 - 7050 A. For the highest redshift targets, 
we needed a redder setting, and used the bluest setting of the 
VPH-Red grism, which covers 5800-9400 A. Our velocity res- 
olution was (Ti nst 67 km s" 1 for both settings. Each night, in 
addition to the science targets, we observed two to four spectro- 
photometric standards covering a range of airmass. Finally, 
over the course of the runs we also observed a large library of 
velocity template stars, consisting mostly of G-M giants. 

The data reduction followed standard procedure^ 
with the exception of pattern-noise removal as dis- 
cussed below. We first performed cosmic-ray 

1 http://iraf.net/irafdocs/spect/ 
2 http://www.ociw.edu/Code/cosmos 



Table 1 . Sample and Observations 



SDSS Name 






Obs. Dale 


t s 


N 




PA 


(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


(7) 


(8) 


SDSS JO15716.92-O05304.6 


0.422 


42.77 


2007-08- 


15 


4500 


1 


180 


88 


SDSS J084135.09+010156.3 


0.111 


AIM 


2007-03- 


18 


7200 




60 


44,74 


SDSS 1110621.96+035747.1 


0.242 


42.59 


2007-03- 


18 


54(10 




60 


60,149 


SDSS J112437.40+045618.8 


0.283 


42.66 


2007-03- 


19 


36(10 


1 




37 


SDSS J 1 14258.39+102747.8 


0.223 


42.81 


2007-03- 


19 


3600 






1 19 


SDSS J122217. 85-000743. 7 


0.173 


42.88 


2007-03- 


18 


7200 




60 


84,174,154 


SDSSJ125334.49-034158.1 


0.239 


42.42 


2007-08- 


14 


1800 




120 





SDSS J135646.10+102609.0 


0.123 


42.77 


2007-03- 


18 


5600 




60 


179,134 


SDSSJ141315.30-014220.9 


0.380 


42.74 


2007-08- 


14 


2700 




600 





SDSS J150719.94+002905.0 


0.182 


42.56 


2007-03- 


18 


4500 




60 


29 


SDSS J205629.85+O05752.6 


0.435 


42.71 


2007-08- 


15 


5400 




600 


109 


SDSS J212653.82+O03553.0 


0.235 


42.27 


2007-08- 


14 


4788 




120 


0.90 


SDSS J221217.92-094407.8 


0.181 


42.26 


2007-08- 


14 


1800 




120 


159 


SDSS J224027.05+0O4347.4 


0.438 


43.02 


2007-08- 


15 


3600 




600 


129 


SDSS J233558.39-0O5O00. 1 


0.439 


42.78 


2007-08- 


15 


3600 




600 


134 



Note. — Col. (1): Official SDSS name. Col. (2): Luminosity distance (Mpc). Col. (3): L|oni] (erg 
s~' ) as measured from the SDSS spectrum of Reyes el al. 2008; used in the object selection. Col. (4): Dale 
of LDSS-3 observation. Col. (5): Total duration of spectroscopic exposure (sec). Col. (6): Number of slil 
positions. Col. (7): Position angle (dcg) ol each slit position, measured E from N. 

removal on each in dividual frame usin g the spectroscopic ver- 
sion of LACosmic ( Ivan Dokkumll200ll) . We performed bias- 
subtraction, fiat-fielding, wavelength calibration, pattern noise 
removal (when necessary; see below) and rectification using 
the Carnegie Observatories reduction package COSMOS0 The 
rectification works extremely well; in the rectified images the 
trace centroid typically migrates by less than two spatial pix- 
els across the full spectrum, with the maximum drift being five 
pixels. COSMOS o ptiona lly performs sky subtraction prior to 
rectification (Kelson 2003), but we have found that for the pur- 
poses of stellar velocity dispersion measurements we achieve 
significantly better sky subtraction using apall in IRAF on the 
rectified, but not sky-subtracted, frames. We have used several 
apertures for spatial extractions, but for most of this paper we 
focus on a 2 "25 spatial extraction, which is well-matched to 
the point-spread function (PSF) for the majority of our o bser- 
vation s. Generally we did not use optimal extraction dHornd 
1 1 986b in order to avoid emission-line clipping. 

Once the spectra were extracted, we performed flux calibra- 
tion, minor shifts to the wavelength calibration and telluric ab- 
sorption c orrections using I DL ro utines, following methods de- 
scribed bv lMatheson et alj (12008). These routines provide con- 
siderable flexibility and control in the selection of bandpasses 
for spectrophotometric calibration used to perform a spline fit to 
the continuum. Small corrections are made, if necessary, to the 
wavelength solutions using a cross-correlation with sky lines, 
and a heliocentric correction is applied. Our spectrophotomet- 
ric standards are typically white dwarfs, and thus the smooth 
continuu m is ideally suited to telluric absorption re moval as 
well (e.g.. lWade & Hornelll988l;lMatheson et alj|2000h . Below 
we show that our flux calibration agrees reasonably well with 
the SDSS by comparing flux measurements of the [O III] line. 

The final spectra have a considerably higher signal-to-noise 
ratio (S/N) than the SDSS spectra. For z < 0.25, the median 
S/N at 5080-5400 A is 150 per pixel, which is typically 4 
times higher than the SDSS, while for z > 0.25, the median S/N 
over the same region is 1 1 per pixel and 8 times higher than the 
SDSS spectra. 

LDSS3, particularly during the March 2007 observation pe- 
riod, had intermittent but substantial pattern noise. Pattern 
noise appears as periodic striping in the data, with a typical 





Rest Wavelength (A) Rest Wavelength (A) 

FIG. 1 . — Example fits to two objects in the sample, spanning the range of redshift and S/N in the sample. Left: Stellar velocity dispersion fits using the best 
weighted fit of an F2, G5, Kl, and K4 star. We plot LDSS3 spectra (thick histogram), best-fit model (red solid), and residuals (thin histogram). Data and models 
are offset for clarity, while regions excluded from the fit are demarcated with grey hatched lines. These fits provide both a measurements of cr* and a continuum- 
subtracted spectrum. Flux density is normalized to the continuum. Right: Fits to the H/3 and [O m] emission lines using multi-component Gaussians; see text for 
details. Total model (thin solid) and model components (dotted) are plotted over the continuum- subtracted data (thick histogram). Residuals are offset for clarity 
and flux density is in units of 10~ 15 erg s _1 cnT 2 A -1 . 



scale of a few pixels and an amplitude (in these data) of 50-100 
counts. The amplitude and phase vary in time and so standard 
calibrations do not remove it but there is a COSMOS task that 
Fourier transforms the data and removes the high-frequency 
noise. Our concern is that the frequency of the pattern noise 
is comparable in scale to that of the expected stellar veloc- 
ity dispersions, potentially leading to biases in our final mea- 
surements. In order to test this possibility, we have performed 
three independent pattern-noise removal runs for the first night, 
which has particularly bad pattern noise. We have measured the 
dispersions independently for each of these spectra. The rms in 
the measured er* values for the four objects range from 3 to 20 
km s" 1 , and, most crucially, we see no systematic differences 
in the measured dispersion for different versions of the pattern- 
noise removal. We conservatively add an additional 20 km s" 1 
uncertainty in quadrature to the measured dispersions that suf- 
fer from pattern noise. 

3. STELLAR VELOCITY DISPERSIONS 

One of the primary uncertainties in interpreting observations 
of narrow-line quasars comes from their unknown BH masses. 
For the median L[oni]= 2 x 10 42 erg s" 1 of the Magellan sample, 
we estimate that the bolometric luminosity is L\, Q \ ps 2 x 10 45 
erg s , which suggests BH masses of > 10 7 M Q if they emit 
at or below their Eddington limits [where we take the Edding- 
ton luminosity to be 1.26 x 10 38 (M B h/M ) erg s" 1 ]. However, 
given the scatter between [O III] and bolometric luminosity, the 
range of bolometric luminosities in the sample, and their (pre- 
sumed) range of intrinsic Eddington ratios, it would be nice to 



have an independent estimate of the BH masses. At the same 
time, it would be very useful to have a mass scale for the host 
galaxies. Although it is indirect, if we assume that the Mbh - 
relation holds, then we can use measurements of bulge veloc- 
ity dispersions both to provide an estimate of the BH masses in 
the sample and as a proxy for the mass of the host galaxy (e.g., 
iHeckman et ail 2004). Here we describe the measurements 
and uncertainties; the BH masses are presented in §5. 

There are a variety of commonly used techniques for mea- 
suring stellar velocity dispersions in galaxies. For computa- 
tional expediency, most of them rely on Fourie r techniques, 
including cross-cor relation ( Tonrv & Davis! 1979 ). the Fourier 
quotient method dSimkinl [19741: ISareent et all 1 1977b . and the 
Fourier correlation quotient ~(lBenderl ll990). It is also possi- 
ble to directl y compare broadened templates with data in pixel 
space (e.g., |Burbidge et al. 1961). Direct-pixel fitting meth- 
ods, while relatively computationally expe nsive, provide many 
benefits over Fourie r techniques (e.g., iRix & White! 11992 ; 



2002; 



van der MareH Il994t iKelson et al.1 l2000t iBarthetal 
Bernardi et al.l 120031) . including the ability both to fit over re- 
stricted regions and to mask complicati ng emission or absorp- 
tion features. As described in detail in Gr eene & Hoi (|2006), 
even faint emission lines can lead to substantial biases in de- 
rived a. t measurements if not properly accounted for. 

Our galaxy models are built from a linear com- 
bination of stellar spectra. The program object is 
de-redshifted using the SDSS redshift, but inevitably 
small velocity shifts are present between the templates 



GREENE ET AL. 



Table 2. Galaxy Properties 



Object 


Bit 


it 


M r 


k 


A f 


'"[Olll] 


r e 


5 


B/T 


0> 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1 111) 


(11) 


JO 157 




4 


-23.16±1.18 


0.50 


0.08 


20.9 


2.6±9.1 


5.7 


1.0 


160±21 


J0841 




2 


-21.70±0.15 


0. 14 


0. 1 2 


21.6 


1.2±0.1 


2.1 


1.(1 


214±29 


Jl 106 




4 


-22.08±0.20 


0.30 


0.16 


20.7 


0.2±0. 1 


3.9 


| 


187±25 


Jl 124 






-23.82±0.15 


0.3 1 


0. 1 5 


22.4 




4.4 




283±21 


J1142 


r 




-22.86±0.15 


0.21 


0.17 


20.7 




3.7 




123±25 


J1222 


r 


4 


-21.10±0.40 


0.1 1 


0.08 


19.8 


0.3±0.4 


3.0 


(1.1 


187±40 


J1253 


r 


2 


-19.82±0.51 


0.17 


0.06 


21.3 


1.8±0.2 


3.9 


0.001 


108±23 


J 1356 


r 




-23.19±0.15 


0.1 1 


0.08 


20.5 




2.3 




206±36 


J1413 


g 


4 


-21.86±0.17 


0.61 


0.26 


23.1 


0.4±0.1 


5.4 


1.0 


118±20 


J 1507 


r 


2 


-20.95±0.18 


0.17 


0.17 


21.2 


0.9±0.1 


3.2 


1.0 




J2056 


r 


4 


-22.03±0.16 


0.44 


0.23 


21.9 


0.5±0.1 


5.8 


1.0 


229±27 


J2126 


i 


4 


-23.45±0.19 


0.04 


0.10 


21.9 


2.9±0.5 


3.8 


0.7 


172±19 


J2212 


r 




-22.86±0.15 


0.12 


0.12 


23.0 




3.1 




133±11 


J2240 


r 


2 


-23.02±0.15 


0.66 


0.19 


20.3 


0.4±0.2 


5.8 


1.0 


314±62 


J2335 


g 


2 


-22.41±0.16 


0.79 


0.13 


21.3 


0.5±0.1 


5.8 


1.0 


119±53 



Note. — Col. (1): Shortened SDSS name (see Table 1 for full name). Col. (2): Filter used for structural 
measurements. Col. (3): Best-fit Sersic index for the bulge-like component from the GALFIT model. Col. 
(4): Absolute magnitude (mag) of the bulge component as derived from our best-fit GALFIT model, with an 
[O ill] correction (Col. 7), correction for Galactic reddening (Col. 6) and a ^-correction (Col. 5) derived using 
the code of Blanton & Roweis 2007. Cases without sensible multi-component models list total Petrosian 
magnitudes and are not assigned a Sersic index, effective radius, or B/T. Col. (5): /^-correction (mag) used 
to correct the absolute magnitudes, filter shown in Col. (2). Col. (6): Galactic reddening correction (mag) 
in the filter shown in Col. (2). Derived from the dust maps of Schlegel et al. 1998 and extinction law of 
Cardelli et al. 1989. Col. (7): Magnitude (mag) of the [O ill] line used to correct the total r-band (z < 0.38) or 
/-band magnitude. Col. (8): Effective radius (") of the bulge-like component of our best-fit model. Col. (9): 
Scale-factor (kpc arcsec -1 ) to convert r e to a physical size. Col. (10): Bulge-to-total ratio for objects with 
sensible GALFIT decompositions. The B/T for J 1253 is very low because there is also a bar component. If 
included as bulge light, the B/T increases to 0.05. Col. (11): Stellar velocity dispersion (km s" 1 ) within 2'.' 25. 

and the target. Thus, each template is shifted to zero velocity, 
diluted by a constant or power-law component C(A), convolved 
with a Gaussian G(A), and multiplied by a polynomial P(X): 

M(\) = P(\){[T(X)®G(X)]+C(\)} (1) 

The nonlinear Levenberg- Marquardt alg orithm, as imple- 
mented by mpf it in IDL (Markwardt 2009), is used to min- 
imize x 2 m pixel space between t he galaxy and the mo del. Our 
code is based on that presented in iGreene & Hoi d2006l) . except 
that the model spectrum is actually a linear combination of stars 
(typically four in this work) with their relative amplitudes as 
additional free parameters. This version of the code was devel- 
oped to model the high-quality nuclear spectr a from the Palo- 
mar spectroscopic survey of nearby gal axies dHo et al. || 19951) 
and is presented and extensively tested in lHo et all d2009i) . Un- 
less otherwise noted we have used the following four stars as 
our template set: HD 26574 (F2 III), HD 107950 (G5 III), 
HD 18322 (Kl III), and HD 131507 (K4 III). Example fits are 
shown in Figure 1. While the primary purpose of the multi- 
template fitting is to achieve reliable fits, we also extract crude 
stellar population information in §7. 

In order to use the full spectral range of the observations, 
we cannot rely only on our own template stars, whose wave- 
length coverage does not go blueward of ^4300 A. Rather, 
we utilize a large library of high-S/N, high -spectral reso lution 
(~ 26 km s" 1 ) stellar templates from Valdes et alj (120041) . It is 
best to use template stars observed with an identical instrumen- 
tal setup as the targets, so that the measured broadening can be 
ascribed solely to internal kinematics in the galaxy. However, 
the Valdes templates not only allow us to include spectral fea- 
tures such as the G-band A4304 A to the highest redshifts, but 
also provide a wide range of spectral types to help minimize the 
impact of template mismatch. We thus use a simple bootstrap- 
ping procedure to remove the impact of differing instrumental 
resolutions on the measured dispersions. First of all, we mea- 
sure the intrinsic broadening of the 
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FIG. 2. — Comparison of stellar velocity dispersion measurements us- 
ing our internal velocity template stars (<TQ ur ) and the Valdes template stars 
(""Valdes) f° r the low-redshift targets. The template stars, while not identical, 
are matched exactly in spectral type to remove ambiguity caused by template 
mismatch. The measured dispersions from the Valdes templates are corrected 
for resolution assuming that the instrumental resolution of our LDSS3 mea- 
surements is <t=67 km s~' , while the instrumental resolution of the Valdes 
stars is <r=26 km s~'. As demonstrated here, the agreement is excellent, which 
supports our assumed conversion factor. 

LDSS3 template stars using the Valdes library. While we did 
not observe identical stars, we can match the stars in spectral 
type nearly exactly. We expect the relative broadening to be 
^rei = (CTldss3- Valdes) 5 = (67 2 -26 2 ) 05 = 62 km s" 1 . We mea- 
sure a mean broadening over 10 stars of 60±7 km s -1 . Dividing 
the sample into the two runs demonstrates the stability of the 
spectrograph; we find 61 ±8 km s" 1 and 59±7 km s" 1 , respec- 
tively. Here and throughout, when we quote errors on mean 
quantities, we are quoting the standard deviation of the distri- 
bution rather than the error on the mean. 

We then confirm that our conversion works properly by mea- 
suring the velocity dispersions of our targets with the Valdes 
templates and comparing the values to those measured from 
our internal templates, again being careful to match the spec- 
tral types of the template stars and the fitting regions. The re- 
sulting dispersions are shown in Figure 2, where the reported 
Valdes dispersions have been reduced by 60 km s" 1 in quadra- 
ture. It appears that our conversion is quite robust. From now 
on we report results based entirely on the Valdes stars with the 
above correction, conservatively adding 7 km s" 1 in quadrature 
to our error budget. The resulting cr„ measurements are shown 
in Table 2. 

In the interest of uniformity, in all cases we use the spec- 
tral region 4100-540 A, which is access ible for every ob- 
servation. We follow IGreene & Hoi (2006) and mask out the 
Mg lb AA5167, 5172,5183 triplet. While the Mg lb triplet is 
one of the highest equivalent width (EW) features in this spec- 
tral region, the well-known correlation between the JM g/Fel ra- 
tio and a* (e.g JO'Connelllll976llKuntschner et al.ll200ll) leads 
to a s ystematic overestimate of a* based on Mg lb dBarth etaTJ 
2002). 
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FIG. 3. — Input-output simulations to estimate uncertainties in a, resulting 
from the finite S/N and continuum dilution. A<r= (<T ut — 0"i n )/o"m. We inves- 
tigate S/N=10 (solid black, circles), S/N=15 (blue short-dashed, triangles), and 
S/N=30 (red long-dashed, squares). In all cases, two-thirds of the continuum 
is a nonthermal component modeled as a constant. Each setting is run with 
10 randomly generated error arrays. We see that there is no net offset even 
at these low S/N ratios and high dilution values, but uncertainties range from 
~ 10% to 20%. 



3.1. SDSS Spectra 

In order to increase our sample size, we have also considered 
all of the targets in the lReves et al.l d2008l) sample. Since most 
of the spectra do not have the required S/N for reliable velocity 
dispersion measurement, we have limited our attention to ob- 
jects with a median S/N of >15 per pixel between 5080 and 
5400 A, and an EW in the Ca K line greater than 3 A. There are 
547 objects with Lpiii] > 10 42 erg s" 1 ; of those, 111 make our 
cut. The SDSS velocity dispersions are measured in an identi- 
cal fashion as the fits described above, using the same mix of 
Valdes template stars. In this case, we assume an instrumental 
resolution of 71 km s" 1 for the SDSS spectra (see discussion in 
iGreene & Holl2006l) . 

The SDSS-derived Z40111] values, our measured velocity dis- 
persions, and [O III] line widths for this additional sample are 
all presented in Table 3. Two objects are unresolved (er*< 
71 km s" 1 ). There are only two objects from the Magellan 
sample for which we can measure reliable dispersions from 
the SDSS spectra, and the resulting dispersions agree with ours 
within the quoted uncertainties. 

3.2. Uncertainties 

Our quoted uncertainties are a quadrature sum of the formal 
errors from the fit and the 7 km s" 1 uncertainty incurred in con- 
verting to the Valdes system. Our typical uncertainties thus de- 
rived range from 10% for the lower-redshift targets to 40% for 
the most distant and luminous targets. 

As the S/N in the spectrum decreases, our ability to recover 
the true dispersion decreases. Also, an increasing contribution 
to the continuum from sources apart from old stars lowers their 
EW and thus the S/N. To a small degree the relative incurred 
error also depends on the true dispersion. 



SDSS Name 






^10111] 


0"|O1I] 


fyorm] 


(1) 


(2) 


(3) 


(4) 


(5) 




SDSS J001609.85-100036.2 


0.167 


180±42 


190± 12 


177±9 


42.16 


SDSS J003535.09-000010.8 


0.174 


162± 25 


141± 24 


197± 21 


42.54 


SDSS J005054.62-000157.4 


0.135 


335± 50 


322± 32 


238± 15 


42.17 



Note. — Table 3 is available in its entirety via the link to the machine-readable version 
above. Col. (1): Official SDSS name, corresponding to RA and Dec. Col. (2): Redshift as 
measured by the SDSS pipeline. Col. (3): (7+ (km s _1 ) measured from the SDSS spectra. 
Col. (4): FWHM[onrj/2.35 (km s" 1 ). Col. (5): FWHM [0 ii|/2.35 (km s" 1 ). Col. (6): L [0 m] 
(ergs- 1 ). 

We run a suite of simulations designed to quantify each of these 
effects by adding Gaussian random noise and a constant contin- 
uum of various levels to one of the template stars that has been 
broadened to various widths. We create ten mock spectra for 
each S/N value, and then measure cr, for each mock spectrum. 
This simple input-output experiment shows that our error bars 
are reasonable. We expect 10% -20% errors for the targets with 
10 < S/N < 30 and 67% continuum dilution (Fig. 3). 

Systematic uncertainties may be introduced in the gas and 
stellar dispersions due to our use of fixed angular rather than 
physical apertures. In the case of the stars, we are typically ex- 
tracting the spectra at ~ 0.5 r e . Since the velocity dispersion 
profiles of bulge-dom inated galaxies are known to be flat (e.g., 
Jorg ensen et al.|[l995h . in general aperture bias should be small 
compared to the measurement uncertainties. The only potential 
exceptions are the disk-dominated galaxies (J 1106, J 1222, and 
J1253) for which disk contamination may be significant. In the 
case of the gas, line widths can change dramatically as a func- 
tion of aperture (typically getting narrower at large radius). Un- 
fortunately, the dependence of gas dispersion on radius varies 
substantially fro m object to object (see, e.g., iRice et aL 2006; 
IWalsh et al. 2008) and there is no physically motivated aperture 
to choose. In future work we will explore the full variation in 
gas velocity dispersion and luminosity as a function of radius. 

3.3. Emission-line Fits to Nuclear Spectra 

A convenient additional product of our velocity dispersion 
fits is a continuum-subtracted emission-line spectrum for each 
target. Since we are interested in the [O II] line when it falls in 
the bandpass, we extend the model out to the bluest available 
pixels although we do not use those wavelengths in the fitting. 
Examples of the continuum-subtracted spectra are shown in 
Figure 1 , and we note that all spectra have been extracted with a 
2'.'25 aperture. We fit the H/3, [O III], and [O II] lines with multi- 
component fits. Briefly, the [O III] AA4959, 5007 lines are mod- 
eled with identical sets of Gaussian components, where the cen- 
troid shifts are fixed to the laboratory value. The H/3 line is fit 
with a combination of two Gaussians (but see §7 for an alternate 
fit). Our prescription follows that described in detail in Greene 
& Ho (2005a), except that unlike that work in general we re- 
quire more than two components to adequately model the lines. 
Furthermore, we do not place any restrictions on the relative 
widths or strengths of the lines. The [O II] AA3726, 3729 lines 
are modeled with identical sets Gaussian components (usually 
two), and the line ratios in the doublet are fixed to one. This pro- 
cedure is necessary to decouple true velocity dispersion from 
the small velocity separation of the two components of the dou- 
blet. 

The line widths and strengths are shown in Table 4. The 
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[O III] line is known to ha ve an asymmetric prof il e, often with a 
prom inent blue wing fe.g JHeckman etalJl981h . |G reene &Hol 
(2005a) find that using a two-component fit and removing the 
blue wing results in a line width that better tracks <r*. How- 
ever, we find that in many cases a two-component model is not 
an adequat e description of these lines. As shown by many au - 
thors (e.g., iGreene & H^|2005a; Bar m et al J 12001 lHoi r2009). 
the FWHM of the line is relatively free of bias from the veloc- 
ity structure at the base of the line. Therefore, we simply take 
(jg=FWHM/2.35, as measured from our best-fit emission-line 
model. 

Although in the bulk of this paper we focus on the extrac- 
tions with a 2. "25 aperture, we have made a variety of other 
extractions, with widths of 3."78, 5."67, and 7"56. The last 
has an effective area comparable to the SDSS fiber (7 arcsec 2 ). 
In Figure 4 we compare the L[omj measurements for our ob- 
jects extracted (typically) within a 7."56 aperture with the SDSS 
measurements. We find overall reasonable agreement, with 
(logL [0 ni],sDSS-logi[oni],our) = 0.01 ±0.2. One object is a sig- 
nificant outlier; our measurement of J2212 is ~ 6 times fainter 
than that derived from the SDSS spectrum. This problem per- 
sists regardless of the flux calibrator used. It is curious, since 
we do not see a similar problem for any other target observed 
on this night. Since the discrepancy is seen in both the con- 
tinuum level and the L[om] value, and since J2212 is a highly 
disturbed and spatially extended system (§4), we suspect that 
we were just unlucky in our slit placement and missed some 
luminous region of the galaxy. None of the conclusions pre- 
sented here are impacted if we use the SDSS luminosity rather 
than ours. If we remove J2212, then the scatter between the 
two sets of measurements drops to 0.15 dex. We note that in 
our own observations L[oni] can vary by as much as 80% be- 
tween different aperture positions, and so we are satisfied with 
this level of agreement. 

4. GALAXY STRUCTURAL MEASUREMENTS 

4.1. Photometric Zeropoints 

We have also reduced the acquisition images (Table 1) us- 
ing standard bias-subtraction and flat-fielding routines within 
IRAF. We do not attempt to remove the pattern noise in the im- 
ages. Each frame is tied to the SDSS photometric system using 
three to six field stars. We use the IRAF package phot to per- 
form aperture photometry on these relatively uncrowded fields, 
and then calculate a weighted average zeropoint based on the 
reported photometric errors in the SDSS photometry. The for- 
mal uncertainties in the final zeropoints are < 0.05 mag. 

4.2. Photometric Decomposition 

The ground-based imaging reaches a typical surface bright- 
ness limit of 24 - 25 mag arcsec" 2 , and it is worth using the 
images to derive bulge-to-disk decompositions for the sample. 
Our goal is to isolate the bulge-l ike component of the host 
galaxy in order to compare with the Faber & Jackson! d 19761) re- 
lation of inactive bulge-dominated galaxies. The Faber-Jackson 
relation can provide a valuable diagnostic as to whether the 
galaxy has reached the relaxed state of local elliptical galaxies. 

We use the t wo-dimensional s urface brightness fitting pro- 
gram GALFIT dPeng et al]|2002f) to model each galaxy. Our 
band of choice is Sloan r, since we have it for the ma- 
jority of the targets. However, for various reasons we are 
forced in a few cases to model either the g- or /-band image. 
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FIG. 4. — Comparison of L[oni] as measured from the SDSS spectra and our 
long-slit data with an aperture (typically) of 7"56. The agreement is reason- 



able, with (L[oui],sdss ~ L \ 



) = 0.01 ±0.2, and the scatter drops to 0.15 



[OIII],our 

dex when the outlier (J2212) is excluded. 

In these cases we convert to r band using the measured color 
from the SDSS Petrosian magnitudes. Throughout, magnitudes 
are quoted in the AB system. 

In comparison with broad-line AGNs the modeling is rela- 
tively straightforward in this case since there is no nuclear point 
source. As a ru le, however, parameter coupling is a major com- 
plication (e.g., Kormendy & Djorgovskil ll989). and thus our 
operating principle is to model each galaxy with the minimum 
number of components. We are particularly interested in the 
bulge components of these galaxies, and we only introduce a 
disk component if visible in the images. In distinguishing be- 
tween bulges and bars for disk-dominated galaxies we use the 
reduced \ 2 to determine the preferred model, but we generally 
cannot constrain more than two galaxy components robustly. 
For the purposes of this paper we have not placed rigorous lim- 
its on the presence of low surface-brightn ess componen ts (cf. 
IGreene et aT1l2008l) . Our basic model is the ISersicI (1 1 968l) func- 
tion: 



I(r) = I e exp 



-b n \ - 

r. 



l/n 



(2) 



where r e is the effective (half-light) radius, I e is the intensity at 
r e , n is the Sersic index, and b„ is chosen such that 



I(r)2Trrdr = 2 1 I(r)2Trrdr. 



We adopt the analy tic approximation for b„ from 
MacArthur et al.1 d2004 . as adapted from ICiotti & Bertinl 
(1999): 



/o 
the 



(3) 



„ 1 4 

■ 2n-- + 

3 405« 



46 



131 



2194697 



25515« 2 1148175« 3 



306907 17750n 4 ' 
(4) 

T he Sersic model reduce s to an exponential profile for n = 1 and 
a Ide VaucouleursI (|l948) profile for n = 4. Bars may be mod- 
eled as ellipsoids with very low a xial ratios. We follow de Jong 
d 19961 ; see also lFreeman| [l966) and model the intensity dis- 
tribution in the bar as a Gaussian. In total for a given Sersic 
component, the model parameters include the two-dimensional 
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centroid, the total magnitude, the Sersic index, the effective ra- 
dius, the position angle, and the ellipticity (all of which are 
constants with radius for a given component in this version of 
GALFIT). The sky is modeled as a constant pedestal offset de- 
termined using multiple boxes placed in empty regions of the 
image, and Levenberg-Marquardt minimization of \ 2 m pixel 
space is used for parameter estimation. 

Even without a nuclear point source, proper modeling of the 
observed light distribution r equire s convolution with a PSF. We 
follow Ravindrana th et al.l ((2006) and use Daophot as imple- 
mented within IRAF to build a PSF model, using no fewer than 
five isolated field stars on the same chip as the target galaxy. 
The PSF thus constructed is both noise-free and directly based 
on the PSF from the current image. 

Two-dimensional modeling provides a few advantages worth 
noting. For one thing, we gain the freedom to assign different 
position angle s and (if necessary) ce ntroids to different com- 
ponents (e.g., Wadadek ar et al.| [l999). Another major benefit 
of GALFIT for these galaxies, which often have nearby com- 
panions, is that we can simultaneously model the surrounding 
galaxies rather than mask them, which provides a considerably 
better model of the background. All fits are shown in Figure 5. 

4.3. Uncertainties 

Uncertainties in the measured parameters are driven primar- 
ily by the PSF model and the sky level, as well as inadequacies 
in the assumed model for the surface brightness profile. In or- 
der to estimate the former, we adopt a single star from the field 
as an additional PSF model. We further run a pair of models in 
which the sky level is increased and decreased by 0.5%, which 
is the typical dispersion in sky values across the field. The final 
error bars represent the dispersion in each parameter from these 
four a lternate runs. Our methodology is described in more de- 
tail in lGreene etai1(l2008l) . 

The hardest uncertainties to estimate arise when our assumed 
model is incorrect. Indeed, in J 1356 and J2212 (see Table 2 and 
Fig. 5), the isophotes are so disturbed due to ongoing merging 
activity that we do not achieve a sensible fit from a simple com- 
bination of one or two Sersic components. In these two cases, 
due to the presence of tidal features and spatially overlapping 
nuclei, our models are not particularly meaningful. Instead, we 
adopt a nonparametric Petro sianl (ll 976b magnitude for clearly 
disturbed systems, as noted in Table 2. Furthermore, although 
the direct light from the active nucleus is extincted by many 
orders of magnitude, some of it escapes along unobscured di- 
rections, scatters off of the interstellar matter in the host galaxy 
and reaches the observer. Both our photometric modeling (see 
below) and spectroscopic fitting (§7) may be affected by this 
component. Finally, we do not have usable Magellan images 
for Jl 124 or Jl 142, and so we adopt SDSS magnitudes in Table 
2. 

As a sanity check, we compare our final galaxy luminosi- 
ties to model magnitud es from SDSS, as well a s to magnitudes 
from the HST survey of Zak amska et al.l (120061) for overlapping 
galaxies. The total magnitudes we measure agree very well 
with the SDSS model magnitudes, with an average difference 
(ours-SDSS) of 0.005 ±0.15 mag, excluding J0157 because 
the photometry is exceedingly uncertain due to the proximity 
of a bright star. We are thus confident that our photometry is 
reasonable, but we adopt an additional 0.15 mag uncertainty in 
the bulge luminosities to account for the systematic differences 
between the two sets of measurements. 



There are two objects in common between this paper and the 
HST sample, Jl 106 and J1413. Again, the total magnitudes 
agree reasonably well, within 0.07 mag in the case of Jl 106 
and 0.4 mag in the case of J1413. (Color differences are derived 
from the LDSS3 spectra, and so strictly speaking only apply in 
the inner regions of the galaxy, and we account only for galactic 
extinction.) While we derive similar B/T ratios for Jl 106 from 
the two images, the derived bulge effective radii are different; 
we find r e = 0.8 ± 0.4 kpc, while it is found to be 3.8 kpc in the 
nearest HST band (F550M). The explanation for this difference 
can be seen in the radial profiles derived from the LDSS3 and 
HST images respectively (Fig. 5). A scattered-light component 
is detected in HST images as a biconical structure with a maxi- 
mum extent of 0"9 and appears as a central brightness peak in 
the one dimensional brightness profile. This component is not 
resolved with the LDSS3 images but may bias our derived bulge 
size to low values. At the same time, the strong dust lanes may 
mimic the transition between bulge and disk in the HST anal- 
ysis, causing the HST bulge measurement to be too high. This 
one example emphasizes that we are potentially biased in our 
bulge sizes and luminosities due to scattered light. We do not 
detect the scattered light signal strongly in the spectrum, pre- 
sumably because of internal reddening in the host galaxy. 

4.4. Galaxy Morphology 

In addition to reliable galaxy luminosities, we also have de- 
rived bulge-to-total ratios (B/T; Table 2) and general morpho- 
logical information for this sample. There are four disk galax- 
ies (Jl 106, J1222, J1253, and J2126) while the rest consist 
of a single bulge-like component, or in the case of J 1356 and 
J2212, train-wrecks. Based on SDSS imaging, Jl 124 appears 
to be extended, but not obviously disky, while J 1142 is com- 
pact. It is interesting to note, furthermore, that a large frac- 
tion of the sample have companions (eight out of 15, most con- 
firmed with the long-slit spectroscopy). Four objects are sig- 
nificantly tidally disturbed, and thus obviously currently inter- 
acting (J0841, J1356, J1222, and J2212). There is no correla- 
tion between AGN luminosity or stellar velocity dispersion and 
ga laxy morphology . Even if we include the HST observations 
of IZakamska et all d2006l) . the median L[oni] for objects with 
and without a disk component are identical to within 0.1 dex. 
On th e other hand, the sample is quite small. IZakamska et al.l 
(2008J) find a hint that L[omj and morphology are correlated, 
in the sense that B/T decreases with AGN luminosity, but for 
objects extending to somewhat lower luminosities. 

4.5. k-corrections and Stellar Populations 

Our ultimate goal is to compare the structural properties of 
the narrow-line quasars to those of inactive galaxies, which in- 
volves translating our observations to match the exact observing 
procedures and evolutionary state of the local samples. Since 
our galaxies span a redshift range of 0.15 < z < 0.45, we must 
assume a model of the intrinsic spectral shape in order to cal- 
culate an effective rest-frame luminosity. We use the publicly 
available program k-corr ect based on the SDSS photometry 
(Blanton & Roweis 2007). It would be preferable to use the 
bulge component alone to derive the ^-corrections, but given 
our heterogeneous coverage in different filters for different sub- 
sets of objects we use the SDSS photometry in the interest of 
uniformity. 

Additionally, we must estimate the contamination from emis- 
sion lines. These objects can have high EWs in [O III], mean- 
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FIG. 5. — Fits are shown for each galaxy, with the exception of Jl 124 and Jl 142, for which the PSF is too poor to use the images. Scale bar indicates 10", 
red arrow points North, and black arrow points to the main target. We show the LDSS3 image (left), the best-fit GALFIT model (middle) and the two-dimensional 
residuals (right), and then the one-dimensional radial profile. The points are the data, the red lines are the bulge components, the blue solid line is the disk compo- 
nent, the magenta solid line shows a bar component, and magenta dotted line is the PSF. The green lines are additional n — 2 Sersic components in J1356, which 
is too disturbed to be well fit by azimuthally symmetric models such as these. The one-dimensional residuals are shown at the bottom. The magenta dotted line 
in J2335 represents an unresolved nuclear component. Note that J0157 is in the wings of a very bright star, making the fitting very uncertain. The sharp features 
present in the models of J0157 and J1507 are due to the artificial truncation of the PSF model. Finally, in Jl 106, we compare the ellipse profiles for the HST-ACS 
FR647M (A c = 7210 A; magenta points) and F550M ( A c = 5581 A; green points) with the LDSS3 r-band image. Our bulge size may be contaminated by a prominent 
unresolved nuclear source presumably due to scattered light. 
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ing that the line emission contributes significantly to the broad- 
band luminosity. From 0.1 < z < 0.4, [O III] contaminates the 
r-band magnitude, while at higher redshift it falls in the i band. 
We calculate the total (extended) [O III] luminosity from our 
widest slit extraction, and convert to an AB magnitude using 
synphot in IRAF. The corrections thus derived, as well as the 
^-corrections for the corrected colors, are shown in Table 2. 

Along the same lines, scattered light from the nucleus may 
contribute a sizable fraction of the luminosity in the blue. We 
have tried to disentangle contributions from both star forma- 
tion and scattered light to the total luminosity (§7). Rerun- 
ning k-correct with the blue continuum removed leads to only 
0.02-0.05 mag of change in the final ^-corrections. Finally, of 
course, stellar populations fade as they age, and the magnitude 
of that fading is dependent on the detailed star formation his- 
tory of the galaxies in question. We will explore this last issue 
in §7. 

5. BLACK HOLES MASSES AND EDDINGTON RATIOS 

It is hard to interpret observations of narrow-line AGNs in the 
absence of direct measurements of BH mass. The targets may 
be 10 8 M BHs radiating close to their Eddington luminosity 
or they may be more massive systems in a low accretion state. 
Their importance to the accretion history of the Universe, the 
interpretation of the observed properties of the host galaxies, 
and comparisons with broad-line objects all depend critically 
on the assumed mass scale. The stellar velocity dispersions and 
inferred range of BH masses thus constitute one of the main 
results of this paper. 

Using the Mbh - c* relation of iTremaine et alj d2002l) . we 
convert into BH mass and in Figure 6 we plot Mbh ver- 
sus L[oni] for the sample. Lines of constant Eddington ratio 
are highlighted, where the bolo metric correction , log L\, \ = 
0.991og L [0 iii] + 3.5 is taken from iLiu et all d2009l) . This cor- 
rection is derived in two steps. First, the conversion between 
Lioni i and Lsmo is calibrated using SDSS broad-line AGNs 
(Rey es et al.l l2008h . and then the relation between L5100 and 
bolometric luminosity is derived based on two recent compi- 
lations of s pectral energy distribu tions (SEDs), iMarconi et al.l 
d2004l) and kichards et al.l d2006l) . In the latter case, we have 
extrapolated the UV-optical slope downward beyond the near- 
infrared minimum in order to estimate the full contribution 
of the big blue bump, but exclude the reradiated mid-infrared 
emission. The estimated scatter is 0.4 dex, accounting for the 
0.2 dex systematic difference between the two chosen bolomet- 
ric corrections, the 0.36 dex scatter between L[om] and contin- 
uum luminosity in broad-line objects, and the 0.05 dex scatter 
in the single-b and bolometric correction from Marconi et al. 
dLiu et alj2009l) . 



As seen in Figure 6, our narrow-line quasars appear to span ~ 
2 dex in BH mass, luminosity, and Eddington ratio. Nominally, 
the Eddington ratios cluster around the value log Lboi/iEdd= 
—0.7 with a dispersion of 0.7 dex (see the best-fit Gaussian fit 
to the distribution in Figure 6b), but it is important to bear in 
mind how the Eddington ratios are obtained. As was mentioned 
above, there is an 0.4 dex uncertainty in the conversion from 
L[om] to Lboi- Furthermore, there is a ~ 0.5 dex uncertainty 
associated with each Mbh measurement which arises from the 
uncertainty in the velocity dispersion. Therefore, the observa- 
tions are consistent with a population of sources radiating close 
to their Eddington luminosities, with the distribution of Edding- 
ton ratios broadened only by observational uncertainty. 

Now let us compare our observations with the masses and 



Eddin gton ratios found for b road-line AGNs at similar red- 
shifts, l^ltoeieretal] d2006l) find a strikingly narrow range in 
Edd ington ratio f or l uminous broad-line AGN s at 0.3 < z < 4. 
She n et alj d2008l) and lGavignaud et aU d2008l) report similar re- 
sults, although the latter paper finds evidence for an increased 
spread in Lbd/^Edd at lower luminosity. It appears that our sam- 
ple is consistent with a similar behavior: rather than spanning 
a wide range in accretion rate, the objects are clustered close to 
their Eddington limits. At face value, not only do the narrow- 
and broad-line obj ects share comp arable space-densities at the 
same luminosity dReves et alj|2008l) . but their distributions of 
mass and Eddington ratio are also surprisingly similar. Prob- 
ably our application of a luminosity threshold is partially re- 
sponsible for the observed narrow distribution in BH mass. The 
Eddington limit effectively determines the lowest-mass BHs to 
be included in the sample, while t he space density of h igh-mass 
BHs declines exponentially (e.g., Hopkins et al. 2009). Never- 
theless, while a variety of correlated errors or selection effects 
may contribute to the apparent narrow distributions observed 
in each population, the point we emphasize here is that the se- 
lections and biases are quite different between the broad and 
narrow-line objects, and thus the similarity in derived distribu- 
tions is surprising. 

There are several steps involved in calculating Mbh and 
Lboi/^Edd, each of which has the potential to introduce sys- 
tematic biases into the measurements. First, there is still 
substantial disagreement about the calculation of Mbh in 
broad-line quasars on the basis of continuum luminosity and 
broad line width (e.g. iBaskin & Laorll2005UCollin et alj |2006b 
iGavignaud et alJl2008t lShen et alj|2008l) . Second, the bolomet- 
ric luminosity correction for narrow-line quasars uses a conver- 
sion from L[oni] to Mb derived for broad-line quasars, and this 
procedure may syst ematically bias the bolometric lumin osities 
towards low values dNetzer et al.ll2006fclReves et al.ll2008l) . Fur- 
thermore, one may question our assumption that the Mbh _ er* 
relation holds for narrow-line quasars. For one thing, they are 
active (i.e., currently a ccreting) BHs and thus may b e offset 
from the local relation dHo et al . 20081 iKim et"aTl l2008). In ad- 
dition, there may or may not be significant redshift evolution 
of the local Mbh _ c* relation that we use to calculate Mbh, and 
both the sign and the magnitude of this chan ge are at present un- 
known (IWoo et alj|2006t IPeng et al.ll2006allbl: iTreu et al.ll2007l: 
lLauer et all 12007b . Each of these effects may potentially bias 
the calculation of Eddington ratios by as much as 0.5 dex. 
The similarity of Eddington ratios of broad- and narrow-line 
quasars, as well as the narrowness of their distribution, is an 
interesting and suggestive result, but we caution that some sys- 
tematic effects may be affecting the comparison. 

6. COMPARISON OF STELLAR AND GASEOUS KINEMATICS 

It has long been known that the [O III] line width 
trace s the bulge velocity dispersion in lo w -luminosity AGNs 
(e.g. iHeckman et all 198 it IWhittldl 1 992alfbl: iNelson & Whittle! 
1996 j). Recently, with the discovery of the Mb h-o* relation 
dGebhardt et al.l2000tlFerrarese & MerrittfeOOOl) . interest in the 
gas velocity dispersion (cr g ) has increased, due to the possibil- 
ity of using it as a proxy for er* in distant or luminous active 
galaxies for which it is prohibitive to measure directly (e.g. 
Nels onll2000l; IWang & Lull200lt iBorosonl l2002t IShields et al 
2003:lGrur>e& Mafliurll2004t lOnken et alj|2004t i Greene & Ho 
2005ai iRice et al.ll200a ISalviander et alj|2007l) . As a result, it 



becomes all the more crucial to understand the physical origin 
of the observed correlation, and to characterize whether it de- 
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FIG. 6. — (a) Sample of narrow-line quasars from Reyes et al. I200JD above a luminosity threshold of t[omj= 10 4i erg s _1 for which we were able to measure 
reliable <r» values. The BH masses are derived indirectly from cr*, assuming the Mbh - ""* relation of Trem aine et al. (2002), while the bolometric luminosities 
are derived indirectly from Z,[onrj. We show the Magellan sample (red open circles), the Gemini sample from lLiu et all 120091 ; blue open squares), and the SDSS 
sample (small filled circles). Double nuclei from the Gemini data are indicated with a cross, (b) The inferred distribution of Eddington ratios for the entire sample, 
with the Magellan data (red open), Gemini data (blue hatched) and best-fit Gaussian (thin solid) overplotted. We note that the uncertainties in this derivation are 
very large; 0.4 dex for the bolometric correction and 0.5 dex in A/bh- Given such large uncertainties, the Lbol/^Edd distribution is strikingly narrow. The central 
values of the Kollmeier et al. (dashed), Gavignaud et al. (long dash), and Shen et al. (dash-dot) distributions are shown for comparison. 



pends on properties of the AGN, such as bolometric luminosity, 
Eddington ratio or radio luminosity. This work increases the 
luminosity range over which such comparisons have been done 
by nearly two orders of magnitude, providing a crucial check of 
the assumption that ct»«ct s . By the same token, given that cr g 
and <t* have now been measured for large samples of weakl y 
active galaxies over a wide range of Hubble types dHoll2009l) . 
it is possible to use the variations in this ratio as a probe of 
the impact of the AGN on the surrounding interstellar medium 
(ISM). 

As shown in Figure 7, there appears to be no correla- 
tion between either FWHM[om] or FWHMjon] and ct» at 
these high luminosities (the Spearman correlation coefficient 
is p = 0.07 with a probability P = 0.4 that no correlation is 
present). We note that while we span a narrow range in lu- 
minosity, we actually span a wide range in <r*, suggesting 
that the lack of correlation is a generic property of luminous 
obscured AGNs, and may well be a generic property of all 
AGNs at such high luminosities. We try to salvage the cor- 
relation by relying on lines with lower ionization potential, 
since many studies have found systematic trends between line 
shap e and the critical density or ionization potential of the lines 
(e.g. iPelat et alJll98lUFilippenko & Halperrill984lFilippenkol 
1985; De Robertis & Osterbrock 1986). Generically speaking 
the high-ionization and/or high-critical density lines originate 
from closer to the BH and are thus typically broader, mak- 
ing low-ionization lines a better choice as proxies to ct* (e.g., 
IGreene & Holl2005iJ: iKomossa &Xull2007l) . Instead, we find 
that the [O II] and [O III] gas widths are strongly correlated 
(p=0.9,P= 10" 5 ; FWHM [ oii]/FWHM [OIII] = 1.0±0.2), but nei- 
ther correlates well with ct* . 

In search of some physical insight into the processes driv- 
ing the gas dispersion, we explore correlations between vari- 
ous AGN properties and <7 g /<7*. For convenience, we define 
the logarithmic difference Aer=log er g -log ct*. Within the 
sample as a whole, (Act) = 0.04 ± 0.2. This is to be com- 



pared with (Act) = -0.15 ± 0. 22 for low-luminosity AGNs in 
early-type galaxies dHoll2009l) a nd (Act) = ± 0-46 f or local 
narrow-line AGNs from SDSS (Greene & Hoi l2005al) . His- 
torically, many papers have found a correlation between ra- 
dio power and TO III] line width (e.g., IWilson & Willis! [l 980; 
IWMttldfl992albt iNelson & Whittldl 19961) . We see no evidence 
for increasing Act or cr g with increasing radio power over a 
range of 10 22 < L R /W Hz -1 < 10 25 at 1.4 GHz. The Nelson 
& Whittle sample does extend an order of magnitude higher in 
radio luminosity, but it is clear that jets alone can not explain the 
sources with ver y large gas dis persions (se e also IGreene & Hoi 
l2005ahlHoll2009b . According to|H3 d2009l) . the deviation of ct s 
from ct* shows a clear dependence on AGN luminosity and Ed- 
dington ratio (see also IGreene & Holl2005al) . in the sense that 
more luminous AGNs have a larger CT g /CT* ratio on average. We 
see no statistically significant correlation between L[omj and 
Act (p = —0.1, P = 0.2). Granted, the dynamic range in lumi- 
nosity is very narrow for this sample alone, but, if anything, 
there is a weak negative trend of Act with luminosit y. Indeed, 
the lower-luminosity data from lGreene & Hoi ((2005a) shown in 
Figure led clearly demonstrate that there is not a well-defined 
correlation between Act and L[om] over a much larger luminos- 
ity range. 

While these results are discouraging in terms of the use of 
FWHM[oni] to approximate ct„, they do tell us in no uncertain 
terms that the ISM in the hosts of these luminous quasars is 
very disturbed. As mentioned above, gas line widths in the cen- 
ters of inactive and weakly active galaxies are well-described 
by vi rial or slightly sub-virial motions due to turbulence dHol 
120091) . In low-luminosity active galaxies, there is some ev- 
idence that A GN activity is, on average, in creasing the gas 
line widths dGreene & Hoi l2005at US [2009), but overall the 
gas still appears to be dominated by virial motion, at least in 
a luminosity-weighted sense. In stark contrast, at the highest 
AGN luminosities, we see no correlation of gas and stellar ve- 
locity dispersions, suggesting that the ISM is completely dis- 
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FIG. 7. — (a) Comparison of the velocity dispersion in warm gas as traced by the [O m] line to that of the stars. Following, e.g, Nelson & Whittle 1 1996) we plot 
the FWHM/2.35 for the gas versus t he stellar velocity dispersion, also measured assuming a Gaussian broadening function (§3). We highlight the Magellan sample 
(open circles) and Gemini sample of lLiu et aH (2009, ; open squares). Error bars in FWHM[oni] are dominated by deviations from the assumption of Gaussianity, 
determined by normalizing the measured full-width at 25% and 75% of maximum, (b) Same as (a) but using [O n] to measure the gas line width, (c) Dependence of 
the deviation of FWHMjonn from <r* as a function of L[om] ■ Formally there is no correlation. Symbols as above, but in grey we have added the lower-luminosity 
sample of Greene & Ho (2005a) to increase the dynamic range. Large filled circles show the median and standard error calculated in 0.5 dex bins for each subsample. 
There is no clear correlation between Act and L[om] even over this large luminosity range, (d) Same as (c), but using FWHM[oii] ■ 
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Table 4. Line Emission and Blue Light 



Object 


log £(ant|,nw 


FWHM |0 []i 


[O nr]/H/3 


[O n]/[0 m] 


loglsiocdct 


log Ls m H$ 


logL SLO o.i lm 


log D„(4000) 


H<5 A 


log Mioo 


(1) 




(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(°) 


(10) 


(11) 


J0157 


42.65 


522 


10.62 


0.12 


42.20 


42.30 


42.7 


1.4 


2.8 


-1.5 


J0841 


42.23 


427 


10.51 


0.027 


42.91 




42.9* 


1.6 


-0.2 




J1106 


42.90 


453 


11.37 


0.060 




42.29 


42.8 


2.7 


-5.5 


-2.5 


11124 


42.12 


335 


7.27 


0.15 


42.19 


42.34 


42.7 


1.5 


3.1 


-1.0 


J1142 


42.59 


567 


10.11 


0.092 




42.67 


42.7 


2.5 




-3.9 


J1222 


42.98 


551 


9.40 


0.035 


43.61 




43.6* 


1.6 


0.3 


-1.1 


J1253 


42.36 


1049 


11.62 


0.15 


43.08 


42.33 


43.1* 


1.4 




-1.6 


J1356 


42.56 




7.19 


0.028 


43.30 


42.40 


43.3* 


2.7 


-5.5 


-1.6 


J1413 


42.36 


469 


7.54 


0.12 




41.61 


42.8 


1.3 


4.1 




J1507 


42.49 


547 


6.79 


0.063 


43.02 


42.78 


42.9 


1.6 


0.8 




J2056 


42.42 


669 


7.97 


QM 




42.45 


42.6 






-1.9 


J2126 


42.14 


440 


9.83 


0.21 






42.5 


1.8 


-0.3 


-2.4 


J2212 


41.47 


428 


8.16 


0.14 


43.15 


41.74 


43.2* 


1.6 


-0.2 


-2.0 


J2240 


43.02 


509 


11.60 


0.24 


42.46 




43.1 


2.7 


-5.2 


-1.8 


J2335 


42.61 


417 


9.28 


0.18 




42.46 


42.9 




3.9 


-1.9 



S-, 



3 



Nole. — Col. (1): Shortened SDSS name (see Table 1 for lull name). Col. (2|: Luminosity (erg s" 1 ) in the [O III] line within ihe 1'.'25 a perl lire. Col. (3): 
Full- width al half maximum (km s _l ) of the [O III] line within Ihe I': 25 aperture, corrected for inslru menial resolution. Col. (4): Ratio of [O [ll| Id Ibi emission 
within Ihe 2"25 aperture. Col. (5): Ratio of [O ll] to [O ill] emission from the SDSS spectra. Col. (6): L5100 (erg s ) when detected in our three -component 
continuum model. Col. (7): Li [mi (erg s _l ) inferred from hroad II ! when delected. Col. (M): Final adopted upper limit 011 contribution to Isiihi (erg s~' ) from 
scattered light, as derived from limits on a scattered hroad 1 1 ' com poll en I. Systems where lire hi lie light is substantially more than can he accounted lor with 
scattered light indicated with : '. These system- may have a substantial star forming population. Col. (9): The D..i4i)()()i index, measured from our best-fit 
continuum model lo ihe SDNS spectra following liakigh el al. t 19991. Col. I III): The I lu\, index, measured from our best iil com in cum model to the SDSS 
spectra follow ing Worlhcy & Oitaviani < I \>'>7 >. Positive values of I h\\ indicate ah-.orplioii. Col. (Ill: The fraction of total galaxy mass farmed ill the last .'(It) 
Myr. based -mi stellar poptilalion ivndtesis. 

7. COINCIDENCE OF STAR FORMATION AND BLACK HOLE 
GROWTH 

Despite decades of research, the prevalence of temporally 
coincident star formation and accretion rema ins one of the 
outsta n ding problems i n AGN phy sics (e.g.. |Bor oson & Ok 



1984t iHeckman et alJ 119951: IC analizo & Stockton 
)03bl ~ 



Kauffmann et al.1 12003bt iHol 12005a ). Generically speaking, 



it seems difficult to inhibit star formation when there is a large 
enough gas supply to fuel a luminous AGN. On the other hand, 
robust measurements of star formation are difficult to obtain, 
particularly with a luminous broad-line AGN spilling light in 
all directions and at all wavelengths. Our hope is to gain new 
insights from the study of obscured but intrinsically quite lumi- 
nous accreting BHs. 

Of course, there are always complications. In our case, al- 
though the primary emission from the AGN is obscured, nev- 
ertheless accretion luminosity manages to contaminate many 
standard star formation indicators. The most reliable instan- 
taneous (i .e. < 10 7 yr) rate diagnostic, Ha luminosity (e.g., 

iKennicutl 1 1 9981) . is excited by accretion power rather than 
young stars in our case. Likewise, blue/UV light is prone 
to contamination from scattered AGN light, whose magni- 



Heckman et al. 


H995tlZakamskaet al.ll2005Ll2006llYoung et al.l 


2009> Liu et al. 


2009). In what follows we show that all of our 



galaxies are bluer than their inactive counterparts, and we use a 
variety of indirect techniques to disentangle the young starlight 
from scattered AGN light. 

We begin by modeling the SDSS colors using stellar popu- 
lation synthesis techniques. As input, we use the SDSS ugriz 
magnitudes, with the [O III] and ^-corrections applied, as de- 
scribed in §4.3. We do not apply any correction for potential 
scattered light, as our goal is to establish whether or not there is 
a blue component apparent in the integrated colors of the galax- 
ie s. The stellar populat ion synthesis code is described in detail 
in IConrov et alJ d2009[). In short, the code utilizes isochrones 
from Padova (Marigo & Girardill2007l) . the BaSeL3.1 stellar li- 
braries with thermally pulsating asymptotic bran ch stars added 
dLeieune et alJll997l [199 8: WesteraetaD[2002), and (in this 
case) a lKroup al (120011) initial mass function. Dust attenuation 
optical depths are allowed to vary separately for ol d and young 
stars, following the model of Chariot & Falll d2000t) . 




1.5 2 
D„(4000) 



FIG. 8. — The D„(4000) and H<5a indices meas ured from our co ntinuum 
fits (see text for details). The locus of points from Kauffmann et al. (2003a) 
is shown in contours, while our measurements are the overplotted red crosses. 
Stellar age increases towards the bottom right of the diagram. The top his- 
togram shows the distribution of D„(4000) from the spectra themselves cor- 
rected for [Ne in], but not corrected for the blue power-law component while 
the points are derived from the spectral modeling. 



The fraction of total stellar mass formed in the last 300 Myr for 
each galaxy is shown in Table 4, with a median value of 1 % at 
a median stellar mass of 10 11 M . Already we see that while 
the BHs are experiencing a major growth episode, the galaxies 
are not. 

The integrated colors show that blue light is prevalent, but 
spectroscopy is require d to determine its origin. For the pur- 
pose of comparing with Kauffma nn et al.l d2003al) . we measure 
the D„(4000) and HSa indices. This two-dimensional diagnos- 
tic is predominantly sensitive to age, and relatively insensitive 
to variations in metallicity or dust. We follow their conventions, 
wherein D„(4000) is defined as t he ratio of fluxes in bands at 
3850-3950 a nd 4000-4100 A dBalogh et alJ fl999l). and US A 
is defined in Wor thev & Ottavianil ( 1997F As a function of 
age, galaxies follow a well-defined sequence in D„(4000)-H(5a 
space, and the locus of SDSS galaxies in this space is shown 
in Figure 8. Unfortunately, because our AGN-powered emis- 
sion lines have such high EWs, we cannot measure a mean- 
ingful HSa index directly from the spectra. We can measure 
D„(4000) after correcting for [Ne III] A3869, and those val- 
ues are shown in the top histogram in Figure 8. [We have 
verified that we can reproduce the values in Kauffmann et al. 
Measuring the 30 most luminous objects in their DR2 catalog, 
we find ([D n (4000) our -D n (4000)Kauff]/D n (4000) K auff) =-0.04± 
0.05.] The median of our sample, uncorrected for any blue 
component, is (D„(4000)) = 1.2. For comparison, the thirty 
most luminous targets in the Kauffmann et al. DR2 sample, 



with a median L[oni] of 10 erg s (~ 0.5 dex less luminous 
than our median luminosity), have (D„(4000)} = 1.3. 

Our first task is to measure the luminosity of the blue con- 
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tinuum. The LDSS3 spectra do not always include the 4000 A 
break, and so we use the SDSS spectra for this measurement. 
We wish to model the various spectral components, but to mit- 
igate degeneracy we restrict the model to nothing more than a 
K star, an A star, and a blue power law whose slope is fixed to 
P\ ex A -1 . In principle the slope of the blue component will de- 
pend on the physical process responsible for it; given the mod- 
eling degeneracies and the uncertainties in internal reddening 
and flux calibration, we do not attempt to measure the true spec- 
tral shape of the blue light. Likewise er» is fixed to the best-fit 
value determined in §3. While the fits are not unique, we be- 
lieve each component is reasonably well constrained: a signif- 
icant A or F star contribution is fixed by higher-order B aimer 
absorption lines, while the EW of the old stellar component 
provides a direct probe of the dilution by blue light. By repeat- 
ing the fitting with a power-law slope of -1 .5, we confirm that 
the derived luminosities are insensitive to the assumed power- 
law slope. More reassuringly, fits performed without an A star 
component yield blue luminosities that agree within 0.05 dex 
of the nominal result. 

These fits provide two useful measurements: one is a blue lu- 
minosity, and the other is a (model-dependent) measurement of 
the D„(4000) and Hi5a indices for the underlying stellar popu- 
lation, as shown by red crosses in Figure 8. The stellar popula- 
tion gets older as one moves to the lower right of this diagram. 
Note that three objects have D„(4000) indices that are larger 
than values seen in real galaxies. Presumably the discrepancy 
arises because in these cases the stellar continuum is modeled 
by a single K star, whereas in real stellar populations there will 
be some contribution from M stars as well as younger popula- 
tions. The D„(4000) and HSa indices derived this way represent 
a lower limit on the young population, while below we will at- 
tempt to evaluate what fraction of the blue light may originate 
in young stars. 

If the blue light is predominantly scattered AGN light, we can 
predict the corresponding reflected broad H/3 component, us- 
ing the well-known tight correlation betwee n optical/UV con- 
tinuu m luminosity and Balmer emission (e.g.JSea rle & Sargent 
1968); we use the recent calibration of iGreene & Hoi d200 5b). 
The expected line fluxes are low (~ 10~ 16 erg s" 1 cm" 2 ), but 
we can at least determine whether or not our spectra are con- 
sistent with a broad compon ent at this level. In the most lu- 
minous narrow-line quasars, iLiu et alj d2009l) robustly detect 
a broad H/3 component in a significant fraction of their sam- 



continuum (see also 


Terlevich et alJl990llHeckman et aljl995l 


Cid Fernandes et al. 


2004). Here, because the lines are much 



weaker, we rather take the indirect approach of comparing the 
expected scattered emission from the continuum fits with a di- 
rect fit to the H/3 line. If the two H/3 fluxes are consistent, then 
we conclude that the blue light may be ascribed exclusively to 
scattered light. 

Since the narrow-line shapes are quite complex and highly 
non-Gaussian, care must be taken to ensure that structure in the 
narrow lines does not masquerade as a bona fide scattered broad 
component. We thus repeat the narrow-line fits as described 
above, except that in this case we fix the narrow components 
of H/3 to that of [O III], and allow the fit to include an addi- 
tional broad H/3 component. We constrain 3000 <FWHMh/j < 
6000 km s" 1 (Liu et al. find H/3fwhm « 4000 km s" 1 ), since oth- 
erwise the broad line will often attempt to fit low-level wiggles 
in the continuum. 



In all cases the [O III] and H/3 line shapes match each other 
extremely well (Fig. 9), and in eight cases we detect a broad 
component (Table 4; Fig. 9 bottom). We convert the rms resid- 
uals into an upper limit on broad H/3 assuming a FWHM of 
4000 km s" 1 . We take as the upper limit on the scattered broad 
component either the detected broad H/3 luminosity or the up- 
per limit thereof, whichever is larger, and then compare with 
the luminosity we expect based on the detected blue contin- 
uum. If the two agree within a factor of 2, then we believe scat- 
tered light can fully account for the blue light, whereas when 
the continuum level is higher, we ascribe the excess to current 
star formation (the five cases are noted with * in Table 4). 

When detected, the blue light accounts for anywhere from 
5% to 60% of the total continuum flux at 4000 A. After the 
blue light is removed, our galaxies occupy a sensible region of 
the H<$a-D„(4000) plane relative to the inactive galaxy popula- 
tion. The median D„(4000) rises from 1.2 to 1.6. Even exclud- 
ing those galaxies with D„(4000) larger than what is seen in 
SDSS galaxies yields an average value of 1.5. In other words, 
while there is clearly an age spread, the corrected index values 
suggest that the galaxy light is dominated by old stellar popu- 
lations in the majority of cases. In five galaxies (J0841, J1222, 
J1253, J1356, and J2212), we find a blue continuum in excess 
of what can be explained by scattered light, which we interpret 
as a star-forming component. It is interesting to note that these 
systems tend to have higher star formation rates inferred from 
their broad-band colors as well. Even more intriguing is the 
observation that all of the highly disturbed systems have an ex- 
cess (star-forming) blue component. On the other hand, we see 
no correlation between the D„(4000) values and host morphol- 
ogy. In summary, these host galaxies are forming < 1 % of their 
stellar mass at the present time. About half of the objects are 
consistent with having only old stellar populations, while the 
others show evidence for ongoing (in the form of excess blue 
light) or recent [low values of D„(4000)] star formation at a low 
level. While we cannot place absolute limits on the current star 
formation rate to any great precision, we argue here, and based 
on the structural parameters investigated below, that the galax- 
ies are effectively fully formed. 

In Figure 10 we show the r-band Faber- Jackson relation 
for the galaxy bulge components, with the most luminous tar- 
gets circled in blue and the disk- or bulge-dominated galaxies 
differentiated as triangles or circles, respectively. Obviously 
there is a large amount of scatter in this diagram, much of it 
measurement-based. It appears that nearly half of the targets 
are offset toward more luminous bulges at a given <r*, but we 
caution that there are substantial systematic biases hidden in the 
bulge luminosities, as discussed below. 

We do not see any correlations between location in the Faber- 
Jackson diagram and AGN luminosity, young stars (as traced by 
the blue light), or intermediate-age stars [D„(4000)]. Partly, the 
sample is just too small, spanning too narrow a range in lumi- 
nosity, to display any such trends. At the same time, a number 
of potential uncertainties may impact the measured luminosi- 
ties. There is firstly the contamination from scattered blue light, 
which boosts the observed r-band magnitudes by as much as 0.5 
mag. Also, the targets span a range of redshifts (0.1 < z < 0.45), 
but we know that over time a collection of stars will become less 
luminous as the most massive stars evolve. Unfortunately, as 
shown by, for example, IConrov et alj d2009l) . the precise mag- 
nitude of this fading is highly uncertain (by ~ 1 magnitude), 
due to our inadequate knowledge of the slope of the initial mass 
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Rest Wavelength (A) Rest Wavelength (A) 

FIG. 9. — (Left) Three example fits used to constrain the blue-light fraction. We fit 3680 — 5400 A with a model including a K star, an A star, and a power law 
(P x oc A" 1 - 5 ), while cr» is fixed to our best fit from above. We plot the SDSS data (solid histogram), the best-fit model (red line), and model components (thin 
black solid line) as well as excluded regions (grey hatched). J1253 and J1507 contain a mixture of all three components, while J1356 requires only a K star and 
a (significant) power-law component. (Right) We show our best fits to the [O m]+H/3 region of our LDSS3 spectra, with narrow H/3 constrained to have the same 
shape as [O m], We show the data (solid histogram), the total model (thin solid line), each model component (dash-dot line), the broad H/3 component (red solid 
line), and the expected scattered broad component (blue dashed line). The latter assumes that all of the blue light is scattered, takes a standard conversion between 
H/3 and L5100 I Greene & Ho 2005b), and assumes an H/3 linewidth of 4000 km s - ' . In J1356 and J1253 there is evidence for recent star formation from the excess 
of blue light compared to the apparent broad component at H/3, while J1507 can be explained purely with scattered light. Flux density in units of 10~ 17 erg s -1 cm -2 
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FIG. 10. — (a)Rest-frame r-band Faber-Jackson relation. We differentiate single n > 2 profiles (open circles), bulge components of disk galaxies (triangles), 
and disturbed systems for which we use the nonparametric Petrosian magnitude (filled circles). Active galaxie s with Z,[om]> 10 42 ' 5 erg s are circled in blue. The 
inactive relation (red solid line), and uncertainties thereof (red dashed lines), are from Desroches et al. ( 2007), and are shown only over the region for which they 
were measured, (b) The same as (a) but with both a correction for the blue light (see text) and an evolution correction of 0.9z mag. 



function around the turnoff mass. Taki ng a value of 0.9z mag 
in r fo r the purpo ses of illustration (e.g.. Ivan Dokkum & Franxl 
1200 U iTreu et al.ll200l IConrov et al.l l2009). and removing the 
detected blue light yields the magenta points, substantially low- 
ering the discrepancy with the local Faber-Jackson relation. We 
are reluctant to provide any additional interpretation, given the 
small sample and substantial uncertainties in the measurements. 

Our basic conclusion, combining the results of the stel- 
lar population modeling, stellar population indices, and struc- 
tural properties, is that for the most part the host galaxies had 
their stellar content in place prior to this epoch of substan- 
tial BH growth. The general result, that we see a range of 
star-forming properties with no clear correlation between host 
galaxy morphology and star formation properties, is in broad 
agreement with previous results for lower-luminosity narrow- 
line AGNs (e.g.. ISchmitt et al.l | l999t iGonzalez Delgado et ail 
l200lMCid Fernandes et al.ll2004l) . 

Most interestingly, the situation appear s rather d i fferen t for 
objects at the highest luminosities. There. IlIu et al. 1 (120091) see 
clear evidence for Wolf-Rayet features in ~ 1/3 of their targets 
(we see none in this sample), which constit ute a clear signa- 
ture of a recent burst o f star formation (e.g..lAllen et al.l 19761; 
lKunth& Sargentlll98ll) . Furthermore. IZakamska et all 00081) 
find very high star formation rates from mid-infrared spec- 
troscopy, for a different sample of luminous narrow-line AGNs, 
rates that are substantially higher than either those seen in in- 
active galaxies of similar mass or those in t ypical broad-line 
AGN s amples of similar luminosity. Now, iKauffmann et al.l 
d2003bl) find compelling evidence that the post-starburst frac- 
tion increases with increasing [O III] luminosity. We may be 
seeing a hint that the galaxy growth rate scales with AGN lu- 



(see also Canalizo & Stockton! 2001 ; 


Vanden Berk et al.ll2006; 


Letawe et al.ll2007t Jahnke et al.ll2007 


). 



HST images show that hosts of luminous narrow-line 
quasars from the SDSS-selected sample are bulge-dominated 



(Zakamskaet ail l2005l) . whereas those of somewhat less lumi- 
nous (by about 0.5 d ex) IR-selected n arrow-line quasars are 
dusty, disky systems dLacv et al.l 120071) . This paper presents 
observations of objects with luminosities very similar to those 
in our HST sample, with similar results: about a third of the ob- 
jects are found in strongly disturbed interacting systems, while 
the rest are found in bulge-dominated galaxies. We previously 
suggested that the difference in morphology may reflect the 
difference in the q uasar fueling mech anisms as a function of 
quasar luminosity (Zakamsk a et al.ll2008l) . but of course this in- 
terpretation is only tentative since the comparison samples were 
selected using very different methods. A detailed study of hosts 
of optically selected narrow-line quasars at lower luminosity is 
needed to understand the origin of the morphological differ- 
ence. 

8. DISCUSSION AND SUMMARY 

We have chosen to explore the host galaxy properties of ob- 
scured active galaxies in particular because of the increased 
contrast between galaxy and nucleus. On the other hand, it is 
our burden to determine whether or not our conclusions can 
be generalized to all accreting BHs at the epoch being probed. 
Rather surprisingly, the distributions of BH mass and Edding- 
ton ratio for our sample of narrow-line AGNs are quite sim- 
ilar to those seen in broad-line samples at similar luminos- 
ity and redshift (e.g.. [kollmeieret all 120061; IShen et all [2008; 
iGavignaud et al.l 120081) . Thi s fact, combined with the obser- 
vation of iReves et all ([2008) that the space densities of the 
broad- and narrow-line objects are comparable at fixed lumi- 
nosity, suggest that the two may originate from similar popula- 
tions. There is another aspect to the story, as illuminated by the 
differences in star-forming properties between the narrow- and 
broad-line objects. 

A variety of recent papers have devised clever techniques for 
disentangling photoionization by stellar and nonstellar sources 
using the line strengths of high- and low-ionization lines, 
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where the former are certainly dominated by the AGN but 
the latter can be si gnificantly cont aminated by star forma- 
tion (e.g..lHoll2005bl:lKim et alJl2006l:IWlar-Martm et alJ2 008: 
Melendez et all l2008h ! These studies all come to the same 
conclusion; star formation is more energetically significant in 
narrow-line than in broad-line AGNs, part icularly at the lu- 
minous end (see also Zakamska et al. 2008). The excess star 
formation in narrow-line objects has been interpreted as evi- 
dence for an evolutionary scenario, in which obscured accre- 
tion and ongoing star formation coexist, and then accretion ter- 
mina tes further star formation , yielding an unobscured obj ect 
(e.gJSanders & Mirabe]||l996tlCanahzo & Stocktonl^OOUlHol 
2005b; Hopki ns et alj2006l) . In this picture some event (usually 
a merger) induces both the star formation and nuclear activity, 
and the broad-line object signifies the final stage in the accre- 
tion cycle. 

Perhaps, on the other hand, there is not necessarily a causal 
connection between the star formation event and the accretion 
event, but when both occur simultaneously the odds are higher 
that the accretio n event will be obscured by star fo rmation in 
the host galaxy dRigbv et al.ll2006t ILacv et alj|2007l) . Our data 
alone clearly cannot distinguish between these two possibilities. 
However, we do see a range of star formation properties across 
even this relatively small sample. Naively, if there were truly 
a causal connection between star formation and nuclear activ- 
ity, one would expect that the majority of our sample would 
show a clear sign of ongoing star formation, even if only in the 
form o f post-starburst sign atures. We look forward to analyzing 
the full R eves et alj d2008l) sample to more definitively quantify 
the star formati on properties as a function of AGN luminosity 
dLiu et alJl2009t) . 

We have seen that luminous obscured active galaxies with 
0. 1 < z < 0.45 and 10 43 < L bo i < 10 46 5 erg s" 1 are composed of 
BHs with masses ~ 10 8 M Q radiating at a high fraction of their 



Eddington limits. We are seeing unambiguous evidence that 
BH accretion impacts the host on a galaxy-wide scale; the host 
galaxy ISM is clearly substantially disturbed relative to less ac- 
tive systems. As a result, we cannot recommend the use of 
gas line widths as a substitute for stellar velocity dispersions at 
such high luminosities, at least not in obscured sources. There 
is a clear blue component in the galaxies, apparent both in the 
overall colors and in nuclear spectroscopy, and in most cases 
the blue light is consistent with being scattered emission from 
the obscured nucleus. Nevertheless, nearly half of the sample 
shows evidence for recent or ongoing star formation, while the 
rest are consistent with predominantly old stellar populations. 
Although star formation activity may rise with vigorous AGN 
activity, there is a wide dispersion at all luminosities, suggest- 
ing that we are far from having a consistent picture of the co- 
evolution of host galaxies and BH growth. 
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